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We present measurements of the branching fraction and CP-vioIating asymmetries in the decay 
B° -^ fo{980)K'^. The results are obtained from a data sample of 123 x lO'* 7(45) -^ BB de- 
cays. From a time-dependent maximum likelihood fit we measure the branching fraction B{B — > 
/o(980)(-^ 7r+7r")A:°) = (6.0 ± 0.9 ± 0.6 ± 1.2) x 10"*^, the mixing-induced CP violation parameter 
S = -1.62t^j^ ± 0.09 ± 0.04 and the direct CP violation parameter C = 0.27 ± 0.36 ± 0.10 ± 0.07, 
where the first errors are statistical, the second systematic and the third due to model uncertainties. 
We measure the /o(980) mass and width to be mj-gfggo) = (980.6 ± 4.1 ± 0.5 ± 4.0) MeV/c^ and 

J-/o(980) 



= (43Ig ± 3 ± 9) MeV/c^ respectively. 



PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



In the Standard Model (SM), CP violation arises 
from a single phase in the three-generation Cabibbo- 
Kobayashi-Maskawa quark- mixing matrix [l|. Possible 
indications of physics beyond the SM may be observed 
in the time-dependent CP asymmetries of B decays dom- 
inated by penguin- type diagrams to states such as 4>K'^ , 
r]'K°, K+K-K°, and /o(980)/i:° 0. Neglecting CKM- 
suppressed amplitudes, these decays carry the same weak 
phase as the decay B" — > J/ipK'^ [3j. As a consequence, 
their mixing-induced CP- violation parameter is expected 
to be -?// X sin2/3 = -77/ x 0.74 ± 0.05 4] in the SM, 
where ry/ is the CP eigenvalue of the final state /, which 
is +1 for /o(980)i4:!?. There is no direct CP violation 
expected in these decays since they are dominated by 
a single amplitude in the SM. Due to the large virtual 
mass scales occurring in the penguin loops, additional 
diagrams with non-SM heavy particles in the loops and 
new CP-violating phases may contribute. Measurements 
of CP violation in these channels and their comparisons 
with the SM expectation are therefore sensitive probes 
for physics beyond the SM. 

In this Letter we present a measurement of the 
branching fraction and CP-violating asymmetries in the 
penguin-dominated decay P° -^ foK^ [5j from a time- 
dependent maximum likelihood analysis. We also mea- 
sure the mass and width of the /o resonance. We restrict 
the analysis to the region of the Tr+Tr"^" Dalitz plot 
that is dominated by the /o and we refer to this as the 
quasi-two-body (Q2B) approach. Effects due to the in- 
terference between the /o and the other resonances in the 
Dalitz plot are taken as systematic uncertainties. 

The data we use in this analysis were accumulated 
with the BABAR detector at the PEP-II asymmetric- 
energy e+e^ storage ring at SLAC. The data sample con- 
sists of an integrated luminosity of 111.2 fb~ collected at 
the T{4:S) resonance ( "on-resonance" ) corresponding to 
(122.6±0.7) x 10^ BB pairs, and 11.8 fb"^ collected about 
40 MeV below the T(45') ("off-resonance"). In Ref. we 
describe the silicon vertex tracker and drift chamber used 
for track and vertex reconstruction, and the Cerenkov de- 
tector (DIRC), the electromagnetic calorimeter (EMC), 
and the instrumented flux return (IFR) used for particle 
identification. 



If we denote by At the difference between the proper 
times of the decay of the fully reconstructed B° -^ foK^ 
(P^Qj.) and the decay of the other meson (P" ), the time- 
dependent decay rate fci^^g is given by 



/Qt..(Ai) 



-\At\/T 

"IT 



l + QtagS'sin(AmdAt) 



QtagCcos(AmdAt) , 



(1) 



where Qtag = 1(^1) when the tagging meson B^^^ is a 
B"{B°), T is the mean P" Ufetime, and Am^ is the P°P" 
oscillation frequency corresponding to the mass differ- 
ence. The parameter S is non-zero if there is mixing- 
induced CP violation while a non-zero value for C would 
indicate direct CP violation. 

We reconstruct P" -^ fo{-^ TT^Tr^)K'^ candidates from 
combinations of two tracks and a K'^ decaying to tt+tt^. 
For the tt+tt^ pair from the /o candidate, we use infor- 
mation from the tracking system, EMC, and DIRC to 
remove tracks consistent with electron, kaon, and proton 
hypotheses. In addition, we require at least one track to 
have a signature in the IFR that is inconsistent with the 
muon hypothesis. The mass of the /o candidate must 
satisfy 0.86 < m(7r+7r~) < l.lOGeV/c^. To reduce com- 
binatorial background from low energy pions, we require 
I cos6'(7r+)| < 0.9, where 6{tt^) is the angle between the 
positive pion in the /o rest frame and the /o flight di- 
rection in the laboratory frame. The K'^ candidate is 
required to have a mass within 10 MeV/c^ of the nomi- 
nal K^ mass [g and a decay vertex separated from the 
P" decay vertex by at least five standard deviations. In 
addition, the cosine of the angle between the K'^ flight 
direction and the vector between the /o and the K'^ ver- 
tices must be greater than 0.99. 

Two kinematic variables are used to discriminate be- 
tween signal-P decays and combinatorial background. 
One variable is the difference AP between the mea- 
sured center-of-mass (CM) energy of the P candidate 
and ■\/s/2, where -^/s is the CM energy. The other 
variable is the beam-energy-substituted mass mES = 
^(s/2 + Pi • p_b)2/P? — Pn, where the P momentum 
Pb and 



p, • PB)yEf - p 

the four- momentum of the initial state (P^, 



Pi) are defined in the laboratory frame. We require 
5.23 < TOES < 5.29GeV/c2 and |A£:| < 0.1 GeV. 

Continuum e^e^ ^^ QQ il — u,d,s,c) events are tlie 
dominant background. To enhance discrimination be- 
tween signal and continuum, we use an neural network 
(NN) to combine four variables: the cosine of the an- 
gle between the B^^^ direction and the beam axis in the 
CM, the cosine of the angle between the thrust axis of 
the -B°e(. candidate and the beam axis, and the zeroth 
and second angular moments Lo,2 of the energy flow 
about the B'^^^ thrust axis. The moments are defined by 
Lj — J2iPi ^ I cos^ip where 9i is the angle with respect 
to the Bjej, thrust axis of the track or neutral cluster i 
and Pi is its momentum. The sum excludes the i3j?gj, can- 
didate. The NN is trained with off-resonance data and 
simulated signal events. The final sample of signal candi- 
dates is selected with a cut on the NN output that retains 
~ 97% (52%) of the signal (continuum). 

The signal efficiency determined from Monte Carlo 
(MC) simulation is (37.2 ± 3.1)%. MC simulation shows 
that 4.7% of the selected signal events are misrecon- 
structed, mostly due to combinatorial background from 
low-momentum tracks used to form the /o candidate. 
In total, 7,556 on-resonance events pass all selection 
criteria. We use MC-simulated events to study the 
background from other B decays. The charmless decay 
modes are grouped into eight classes with similar kine- 
matic and topological properties. The modes that decay 
to the Tr^Tr~ K'^ final state are of particular importance 
since they have signal-like AE and to.es distributions 
and their decay amplitudes interfere with the f^K^ decay 
amplitude. Among these modes are p^K^., /o(1370)i^°, 
/2(1270)ii'3, K*^-!!^ (including other kaon resonances 
decaying to iirgTr"*"), and non-resonant 'k^-k~K^ decays. 
The inclusive charmless 'n~^'n~K^g branching fraction 
(23.4 ± 3.3) X 10"^ IIk together with the available ex- 
clusive measurements jj], are used to infer upper limits 
on the branching fractions of these decays. Along with 
selection efficiencies obtained from MC, these branch- 
ing fractions are used to estimate the expected back- 
ground. The charmed decays B^ — > D~tt^ — > K'^tt~tt^ 
and B^ -^ D^-n^ -^ K^tt^ti^ contribute significantly to 
the selected data sample. Each of these modes is treated 
as a separate class. Two additional classes account for 
the remaining neutral and charged 6 — > c decays. In the 
selected data sample we expect 45 ± 15 charmless and 
128 ± 74 6 ^ c events. 

The time difference Ai is obtained from the measured 
distance between the z positions (along the beam direc- 
tion) of the i?°(,(. and B^^^ decay vertices, and the boost 
/37 ~ 0.56 of the e+e^ system [aQ- To determine the 
flavor of the B^ we use the tagging algorithm of Ref. p| . 
This produces four mutually exclusive tagging categories. 
We also retain untagged events in a fifth category to im- 
prove the efficiency of the signal selection. 

We use an unbinned extended maximum likelihood fit 



to extract the foK^g event yield, the CP parameters de- 
fined in Eq. |^, and the /o resonance parameters. The 
likelihood function for the N^ candidates tagged in cate- 
gory k is 



Lk = 



i=l '^ 



^kVh + Nc,kV^ 



i.k 






(2) 

where N^^ is the sum of the signal, continuum and the 
ns -B-background yields tagged in category k, Ns is the 
number of JqK^ signal events in the sample, efe is the 
fraction of signal events tagged in category fc, Nc,k is 
the number of continuum background events that are 
tagged in category fc, and NB,j£j.k is the number of 
B-background events of class j that are tagged in cat- 
egory k. The _B-background event yields are fixed pa- 
rameters, with the exception of the D^tt^ yield. Since 
B^ -^ D'TT^ events have a characteristic distribution 
in cos0(7r~''), well separated from continuum and JqK^ 
events, the D~'k'^ is free to vary in the fit along with the 
signal and continuum yields. The total likelihood C is 
the product of the likelihoods for each tagging category. 

The probabihty density functions (PDFs) V^ , V^ 
and Vff.^ for signal, continuum background and B- 
background class j , respectively, are the products of the 
PDFs of six discriminating variables. The signal PDF 
is thus given by "Pf = 'P'^(toes) • V^{AE) ■ T'f (NN) • 
V^{\ cos6l(7r+)|) • 7'^(TO(7r+7r-)) • Pf (At), where Pf (At) 
contains the time-dependent CP parameters defined in 
Eq. QJ , diluted by the effects of mis-tagging and the At 
resolution. 

The signal PDFs are decomposed into distinct distri- 
butions for correctly reconstructed and mis-reconstructed 
signal events. The fractions of mis-reconstructed signal 
events per tagging category are estimated by MC simu- 
lation. The Toes, AE, NN, | cos6i(7r+)|, and TO(7r+7r") 
PDFs for signal and B background are taken from the 
simulation except for the means of the signal Gaussian 
PDFs for Toes and AE, and the mass and width of the 
/o, which are free to vary in the fit. We use a relativistic 
Breit-Wigner function to parameterize the /g resonance. 
The At-resolution function for signal and i?-background 
events is a sum of three Gaussian distributions, with pa- 
rameters determined by a fit to fully reconstructed B^ 
decays |9|. The continuum At distribution is parame- 
terized as the sum of three Gaussian distributions with 
two distinct means and three distinct widths, which are 
scaled by the At per-event error. For the i?-background 
modes that are CP eigenstates, the parameters C and S 
are fixed to and ± sin 2/3, respectively, depending on 
their CP eigenvalues. For continuum, four tag asymme- 
tries and the five yields Nc,k are free. The signal yield, 
5", C, and the /o mass and width are among the 41 pa- 
rameters that are free to vary in the fit. The majority of 
the free parameters are used to describe the shape of the 



TABLE I: Summary of systematic uncertainties. The uncer- 
tainties on rufg and Tf^ are in units of MeV/c^. 



Error Source 



5" 



C 



B 
xlO" 



"t/o 



fa 



Fitting Procedure 0.06 0.07 0.26 0.5 1.0 

B-background 0.05 0.05 0.30 0.2 3.0 

At Model 0.03 0.02 0.03 0.0 0.1 

Tagging Fraction 0.04 0.02 0.01 0.0 0.2 

Signal Model 0.02 0.01 0.03 0.1 0.2 

DCS Decays 0.01 0.04 0.00 0.0 0.0 

Arrtd and r 0.01 0.00 0.01 0.0 0.1 

Sub-total a09 OTO 040 05 3^ 

Q2B Approximation 0.04 0.07 1.21 4.0 8.5 

continuum background. 

The contributions to the systematic error on the sig- 
nal parameters are summarized in Table Q] To estimate 
the errors due to the fit procedure, we perform fits on a 
large number of MC samples with the proportions of sig- 
nal, continuum and _B-background events measured from 
data. Biases of a few percent observed in these fits are 
due to imperfections in the likelihood model such as ne- 
glected correlations between the discriminating variables 
of the signal and B-background PDFs and are assigned 
as a systematic uncertainty of the fit procedure. The er- 
ror due to the fit procedure includes the statistical error 
on the bias added in quadrature with the observed bias. 
The expected event yields from the _B-background modes 
are varied according to the uncertainties in the measured 
or estimated branching fractions. Since _B-background 
modes may exhibit CP violation, the corresponding pa- 
rameters are varied within their physical ranges. We vary 
the parameters of the At model and tagging fractions in- 
coherently within their errors and assign as a systematic 
error the quadratic sum of the observed changes in our 
measured parameters. The uncertainties due to the sim- 
ulated signal PDFs are obtained from a control sample 
of fully reconstructed S" -^ D~{-^ Xg7r~)7r+ decays. 
The systematic errors due to interference between the 
doubly-Cabibbo-suppressed (DCS) b — > ucd amplitude 
with the Cabibbo-favored b — > cud amplitude for tag-side 
B decays have been estimated from simulation by vary- 
ing freely all relevant strong phases \M^- The errors as- 
sociated with Arud and r are estimated by varying these 
parameters within the errors on the world average [3|. 

The systematic error introduced in the Q2B approx- 
imation by ignoring interference effects between the /o 
and the other resonances in the Dalitz plot (as listed 
earlier) is estimated from simulation by varying freely 
all relative strong phases and taking the largest observed 
change in each parameter as the error. While the system- 
atic errors due to interference are comparable to the sta- 
tistical error for the branching fraction and the /o mass 
and width, they are small compared to the statistical er- 
ror for S and C. 

The maximum likelihood fit results in the foK*^ event 




FIG. 1: Distributions of (clockwise from top left) AE, tties, 
cos6(7r^)|, and m(7r'''7r~) for samples enhanced in foKs sig- 
nal. The solid curve represents a projection of the maximum 
likelihood fit result. The dashed curve represents the contri- 
bution from continuum events, and the dotted line (middle) 
indicates the combined contributions from continuum events 
and B backgrounds. For presentation purposes, the region 
0.765 < I cos0(7r^)| < 0.81 has been removed to suppress the 
contribution from D'tt"*" events. 

yield Ns = 93.6±13.6±6.3, where the first error is statis- 
tical and the second systematic. The branching fraction 
corresponding to this yield is 

B{B° ^ /o(980)ifO) X s{fo{980) ^ ^+71^) = 
(6.0 ±0.9 ±0.6 ±1.2) X 10"^ 



where the first error is statistical, the second systematic, 
and the third accounts for the model dependence in the 
quasi-two-body approximation. The systematic error in- 
cludes an uncertainty of 8.2% from differences between 
data and MC in tracking, particle identification and K'^ 
detection efficiencies. Figure ^ shows distributions of 
AE, TTiEs, |cos0(7r+)| and m(7r+7r~), that are enhanced 
in signal content by cuts on the signal-to-continuum like- 
lihood ratios of the other discriminating variables. For 
the CP-violation parameters, we obtain 

S = -1.62+°;^^ ±0.09 ±0.04, 
C = 0.27 ±0.36 ±0.10 ±0.07. 

The time-dependent distributions and asymmetry 
^s'V-B" ^ {Ngo — N go)/ (N go +N go) i'a the tagged events 
arc represented in Fig. 13 The model-dependent mass and 
width of the /o are found to be 

rrif, = (980.6 ± 4.1 ± 0.5 ± 4.0) MeV/c^, 

T/o = (43lg^±3±9)MeV/c^ 

These results are in agreement with previous mass and 
width measurements [3, llll| . 

In summary, we have presented measurements of 
the branching fraction, resonance parameters, and CP- 
violating asymmetries in B^ — > /o(980)(^ 7r"'"7r~)i^g de- 
cays. Our results for S and C are consistent with the 
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FIG. 2: The signal enhanced time distributions tagged as 
Stag (top) and -Btag (middle), and the asymmetry, Ago /go 
(bottom). The solid curve is a projection of the fit result. 
The dashed line is the distribution for continuum background 
and the dotted line is the total B- and continuum-background 
contribution. 

Standard Model at the l.Tcr and 0.8cr levels, respectively. 
The result for S is 1.2(t from the physical limit, and the 
hypothesis of no mixing-induced CP violation is excluded 
at the 2.7(7 level. 
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